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SUMMARY 

Circular membranes with fixed peripheral edges, subjected to arbitrary 
axisymmetric loading were analyzed. A single governing differential equation 
in terms of radial stress was used. This nonlinear governing equation was 
solved using the finite difference method in conjunction with Newton-Raphson 
method. Three loading cases, namely (a) uniformly loaded membrane, (b) a mem- 
brane with uniform load over an inner portion, and (c) a membrane with ring 
load, were analyzed. Calculated central displacement and the central and edge 
radial stresses for uniformly loaded membrane, agreed extremely well with the 
classical solution. 
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INTRODUCTION 


Composite laminates are being extensively used in aerospace applications 
because of their high strength to weight ratios. Because of their widespread 
use, these laminates in addition to several types of loads are subjected to 
impact loads. Due to these impact loads composite laminates suffer visible 
and invisible damage. Damage characterization and residual strength evalua- 
tion of these impact-damaged composite laminates is of utmost importance. 

During low-velocity impact on thin circular 8-ply quasi-isotropic lami- 
nates, the laminates exhibit large deflection behavior. The displacement in 
the laminates can exceed up to one laminate thickness [1]. Because of these 

i 

large deflections, the laminates exhibit mid-plane stretching. The mid-plane 
stretching contributes to the membrane action. Therefore, a damaged composite 
laminate can be analyzed by decomposing the problem into two parts. The first 
part consists of a thin plate with shear and flexural stiffness but no mid- 
plane extenslonal stiffness. The second part consists of a membrane with mid- 
plane extensional stiffness but no shear and flexural stiffness. The results 
of these two types of analyses will yield complete knowledge of displacements 
and stresses in a quasi-isotropic laminate due to low-velocity impact. This 
knowledge is needed for the damage characterization and residual strength 
evaluations of these laminates. 

The first part of the analysis mentioned above is the well known small 
deflection analysis of isotropic plates [2]. The second part, which is the 
analysis of isotropic membranes with arbitrary axisymmetric loading, however, 
needs further attention. The problem of an isotropic membrane subjected to 
surface and edge loads has been studied by many investigators [3-10] . Hencky 
[3], Dickey [4], and Shaw and Perrone [5] determined the deflection of a 
uniformly loaded membrane. Kao and Perrone [6] used a nonlinear relaxation 
method to obtain deflections in a uniformly loaded membrane. Goldberg and 
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Plfko [7] and Weinitschke [8,9] employed power series approaches to obtain the 
solutions for annular membranes. In addition to the power series method 
Weinltschke [8,9] presented an Integral equation approach to the solution of 
annular membranes subjected to surface and edge loads. Callegari and Reiss 
[10] studied the axisymmetric deformations of a circular membrane subjected to 
arbitrary normal pressure by using the shooting method. The reason for the 
limited numerical solutions for arbitrary loaded membranes may be due to dif- 
ficulties in using the conventional approaches used by earlier investigators. 
One of the difficulties may be due to the large differences in magnitudes of 
the radial and transverse displacements and the satisfaction of the two 
governing equations involving these displacements. 

The purpose of this paper, therefore, is to present an alternate formula- 
tion for the analysis of circular isotropic membranes subjected to arbitrary 
axisymmetric loading. In contrast to the earlier governing equations involv- 
ing radial and transverse displacements, in this formulation a single non- 
linear governing differential equation is used in terms of the radial stress. 
The nonlinear differential equation was then replaced by a set of nonlinear 
algebraic equations using difference quotients. Then by using the Newton- 
Raphson method, these nonlinear algebraic equations were numerically solved 
to obtain stresses and displacements. This type of formulation involving a 
single governing equation avoids difficulties encountered with the differing 
magnitudes of radial and transverse displacements. 

This approach is used to analyze (1) a uniformly loaded circular 
membrane, (2) a membrane with uniformly distributed load over the inner 
portion, and (3) a ring loaded circular membrane. The results for the radial 
and tangential stresses as well as the radial and transverse displacements are 
presented for each case. The results for the uniformly loaded circular mem- 
brane are compared with those from the literature. 
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LIST OF SYMBOLS 


a radius of membrane 

E Young's modulus 

h thickness of membrane 

m number of regions 

n node number 

N r radial tension force 

Ng tangential tension force 

p(r) pressure at any radius r 

r,0,z cylindrical coordinate system 

Ac elemental length 

u radial displacement 

u normalized radial displacement 

w transverse displacement 

w normalized transverse displacement 

e_, e Q radial and tangential strain 
r y 

o , a. radial and tangential stress 

r u 

o_ , a Q normalized radial and tangential stress 
r u 

v Poisson's ratio 

ANALYSIS 

Membrane Configuration and Strain-Displacement Relations 
Figure 1 shows an axisymmetrically loaded circular membrane with fixed 
peripheral edge, with thickness h and radius a. Since the deflection 
surface is axisymmetrical, the displacement can be resolved into two com- 
ponents: (1) a component u in the radial direction, and (2) a component w 
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perpendicular to the plane of the membrane* From large deflection theory [2] , 
strain-displacement relations for an isotropic membrane are 



(la) 


and the corresponding stresses are, 



E 

(1 - V 2 ) 


E 

(1- v 2 ) 



X 


(lb) 


Note that the assumption of large deflection but with small strains is 


made in this formulation. 


Equilibrium Equations 

The equation of equilibrium in the radial direction is 
dN 

N - N. + r -j-=- « 0 (2) 

r 0 dr 

If N_ and N 0 in Eq. (2) are replaced by oh and o Q h, 
respectively, and the thickness of the membrane is assumed constant, the 
equation of equilibrium in radial direction is, 

do 

a r ■ °e + r sr - 0 < 3 > 


The equation of equilibrium in the direction perpendicular to the plane 
of the membrane, with general axisymmetric loading p(r) can be written as 

r 


2in:ho 4— + / p(g) 2*5 d£ » 0 

r ar o 


( 4 ) 
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The stress and strain displacement relations (Eq. (1)) and equilibrium 

equations (Eqs. (3) and (4)) when combined form four nonlinear partial 

differential equations with four unknowns a , o a , u, and w. 

r 0 

Derivation of a Governing Equation 

Using the stress and strain displacement relations (Eq. (1)), the radial 
displacement u was expressed as 


u "f (®e " vc, r) 


(5) 


The term ^ was obtained by differentiating the radial displacement u, 
with respect to r. 


du l f 
dr " E 




( 6 ) 


The strain in the radial direction was expressed in terms of the 

radial and tangential stresses and by using the strain-displacement 
relations (Eq. (1)), 



By using ^ from Eq. (6) in Eq. (7) the following relationship was 
obtained: 



From the equation of equilibrium in radial direction (Eq. (3)), 
written as 


( 8 ) 



was 


da 


dr 



( 9 ) 


Substituting Eq. (9) in Eq. (8) and rearranging the terms the following 
relationship was obtained: 
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Equation (14) is a nonlinear differential equation in terms of the radial 
stress that governs the large deflection response of the membrane with arbi- 
trary axisymmetric loading. Several investigators obtained the governing 
equation in terms of the radial stress a^,, similar to Eq. (14). Dickey [4] 
and Weinitschke [8,9] presented the differential equation for uniformly dis- 
tributed load. Callegari and Reiss [10] obtained the differential equation 
for a membrane with an arbitrary axisymmetric loading. 
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SOLUTION METHOD 


To solve the governing equation (Eq. (14)) for the radial stress, Dickey 
[4] used integral equation method, Weinitschke [8,9] used integral equation 
and power series approaches and Callegari and Reiss [10] used the shooting 
method. In contrast, here a numerical method of solution is proposed. The 
nonlinear differential equation was replaced by a set of nonlinear algebraic 
equations using difference quotients. Then using the Newton-Raphson method 
[11], these nonlinear algebraic equations were solved numerically to obtain 
the stresses and displacements. The details of the solution method are as 
follows* 

The solution domain was discretized into m regions and (m + 1) nodes. 
Denoting a^.(n) as the radial stress at the nth node, the governing equation 
(Eq.. (14)) was rewritten as, at any node n: 



(15) 


To simplify the evaluation of the integral in the above equation, the 
following assumption will be made. Consider an ith region with 1 < r < r^ 
Although the applied pressure varies within the region r^ < r < r^ the 
pressure will be assumed to be uniform in this region with a value p ± . The 
magnitude of p ± is assumed equal to the value of the pressure at the mid- 
point of this region, i.e., at r ™ (r^ + r^_^)/2. As the number of regions 
in the model become large, the size of each region reduces and hence the 
variation of the pressure within each region also becomes insignificant. With 
this assumption the governing equation (Eq. (15)), reduces to 
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First and second derivatives of the radial stress from the governing 
equation (Eq. (16)) were replaced by 

do. cr r (n + 1) - o r (n - 1) 

dr nth node 2 Ar ~~ 

and 

d 2 cr. a r (n + 1) - 2a r (n) + cj.(n - 1) 

dr 2 nth node (Ar) 2 

For the circular isotropic membrane with axisymmetric loading the 
boundary conditions were: 

(1) Both the radial displacement u and the transverse displacement w 
equal zero at the fixed edge (r * a). 

(2) The radial displacement u and the slope equal zero at the 

center (r ■ 0) . 

Since the governing equation (Eq. (14)) was derived in terms of the 
radial stress, the transformation of boundary conditions was done by using 
stress and strain-displacement relations (Eq. (1)), and equilibrium equations 
(Eqs. (3) and (4)). When u «• 0 and r ■ a are substituted into Eq. (5) for 
radial displacement, the boundary condition (1) above can be transformed to 

Oq ■ vo r , at r » a (17) 
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Using this boundary condition In the equilibrium equation In radial direction 
(Eq. (3)), (Eq. (17)), was further transformed to 

da r 

0^(1 - v) + r °0, at r - a (18) 


The boundary condition (2) above is a statement of symmetry about r » 0, the 
center of the membrane* This symmetry condition can be expressed in terms of 
the radial stresses, 
do r 

”0, at r » 0 (19) 

da r 

The derivative was replaced by the finite difference quotient as 

4 

4cr r (2) - 3cr r (l) - (7.(3) - 0, at r » 0. 


Using the governing equation (Eq. (16)), 


dcr r (n) 

dr 


+ r 


d 2 cr r (n) 


n 


dr J 



iMiM 

r n h o r (n) 


and boundary conditions (Eqs. (18) and (19)) in the form of finite difference 

quotients at each node, m + 1, nonlinear algebraic equations were obtained. 

These algebraic equations contained m + 1 unknowns, viz., a r (l) ... a (m + 1) . 

By using the Newton-Raphson technique [11] , these (m + 1) equations were 

solved to obtain the radial stress at each node. 

The equation of equilibrium in the radial direction (Eq. (3)), was used 

to obtain the tangential stress, <jg, at each node. The radial 

dw 

displacement, u, and the slope, -j^-, at any node were obtained by using 
Eqs. (5) and (12), respectively, and the a and o a values at that node. 

Then using the boundary conditions, 
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w ■ 0 at r - a. 


dw 

dr 


0 at r 


0 


and the finite difference representation of the slope, 4^, at all the 

dr 

(m + 1) nodes, the transverse displacements, w, were determined* 


RESULTS AND DISCUSSION 

In this section, first a convergence study for the above outlined method 
is presented. Then, the present method is illustrated for circular membranes 
with different loadings, shown in figure 2: 

(1) Uniformly loaded circular membrane. 

(2) A membrane with uniformly distributed load over the inner portion. 

(3) A membrane with a ring load. 

For all the cases analyzed a Poisson's ratio of 0.3 was assumed. The 
stresses and displacements were expressed in the dimensionless forms, using 
Hencky's normalized formulation as. 



w 


w 






a e 


(&) 


171 


li 



Convergence Study 


To study the convergence of the present method, the circular membrane was 
Idealized into m number of regions with (m + 1) nodes, (where nodes are 
numbered from the center to the outside). The number of regions, m, used in 
this convergence study were 8, 16, 32, 60, and 64. Figures 3(a) and 3(b) 
present the relative errors in the normalized maximum deflections and 
stresses, for a uniformly loaded membrane and for a membrane loaded uniformly 
over the region 0 < — < 0.5. The solution shows rapid convergence and about 

£L 

60 regions were found to be necessary for a membrane loaded uniformly over the 
r 

region 0 < — < 0.5, to yield a solution which is within 0.001 percent of the 

3 . 

converged solution. In contrast much fewer than 60 regions were necessary for 
the uniformly loaded membrane. However, a 60 region idealization is used in 
the analysis and all the results are presented for this idealization. 

The Newton-Raphson method [11] used here needs initial o r values for 
the iterative process. The sensitivity of the present method to these initial 
values was studied. Table 1 presents the number of iterations needed to 
achieve convergence for a range of initial "cr r values for the three configu- 
rations studied. The initial "cr r values were as low as three orders of mag- 
nitude smaller than the maximum converged value of and as high as three 

orders of magnitude higher. For the uniformly loaded membrane with each of 
these initial values the present method converged to the same solution. The 
other two configurations showed similar convergence, thus showing insensi- 
tivity to the initial o^. values. When the initial a r values were farther 
away from the converged values the number of iterations needed were around 40 
compared to about 8 when the initial values were closer to the converged 

\ values. 
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Uniformly Loaded Circular Membrane 
The first problem analyzed was that of a uniformly loaded circular 
membrane, for which classical solution [ref* 3] exists* The uniform loading 
was represented by setting the magnitude of loading terms p^ through p 61 
equal to unit values. Using these values in Eq. (16), the governing equation 
for the membrane was solved by using the procedure outline earlier. The 
values of u • and w displacements and stresses cr r and <jg are presented 

in figures 4 and 5, respectively. A comparison of the present solution with 

/ 

Hencky's classical solution [3] and Kao and Perrone [6] nonlinear relaxation 

\ 

method solution is shown in Table 2. The central displacement as well as 
central and edge radial stresses obtained by present method are in excellent 
agreement with the earlier reported results as shown in Table 2. 

A Membrane With Uniformly Distributed Load Over the Inner Portion 

The second example is that of a membrane with uniformly distributed load 

over an inner portion. A particular case of loading over the region 

® ^ ^ was analyzed. In the corresponding solution, this loading was 

represented by setting the magnitude of loading terms p^ through p 6 equal 

to unit values and terms py through p^^ equal to zero. Again using the 

governing equation (Eq. (16)) with the loading terms discussed before, the 

solution was obtained for stresses and displacements at each node. The 

corresponding values of u and w displacements and a and a n stresses 

r 0 

are shown in figures 6 and 7, respectively. 

A Membrane With a Ring Load: 

The last problem considered here is one where the membrane carried a 
uniformly distributed ring load. As a specific case the ring load was assumed 
to be spread over the region 0.5 < — < 0.6. Since loading was considered over 
the region 0.5 < ~ < 0*6, in the corresponding solution, this loading was 
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represented by setting the magnitude of loading terms through p^ equal 

to zero, P32 through P37 equal to unit values; and P33 through pg^ equal 
to zero. Using these values In Eq. ( 16 ), the governing equation for the 
membrane with the ring load was solved by using the procedure outlined earlier. 

The corresponding values of u and w displacements are presented In 

figure 8 and stresses a and o fl are presented in figure 9 . Figure 8 shows 

that the transverse displacement w is constant up to the ring load and then 

starts decreasing and becomes zero at the clamped edge. In contrast, the 

radial displacement u is zero at the center, increases gradually until it 

reaches the peak value In the loading region and then starts decreasing and 

becomes zero at the fixed edge. On the other hand, the normalized stresses 

a and o a are of the Identical magnitude in the unloaded region 
r 0 

0 < — < 0.5 and then decrease for larger values of r. 
a 

As pointed out in the introduction, when u and w differ by large 
amounts a simultaneous method of solution presents difficulties. As the 
present method does not use u and w as parameters, rather uses a single 
parameter, 0^ these difficulties are avoided. Therefore, it is interesting 
to compare the differences in magnitudes of u and w values for various 
problems analyzed. For a membrane with uniformly distributed load u and w 
are of the same magnitude, whereas in the case of membrane loaded over an 
inner portion and for the ring loaded membrane, the displacements differ by 
two and one order of magnitude, respectively. 

Because u and w do not differ by large amounts in the case of uni- 
formly loaded membrane, one would expect the simultaneous solution to be effi- 
cient and feasible. Indeed, it is so as demonstrated by Kao and Perrone [6] . 
Because u and w differ by large amounts for the arbitrarily loaded mem- 
branes, one would expect difficulties with simultaneous solution method. This 
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may be the reason for the limited numerical of the solutions for these two 
loading cases. The present method on the other hand avoided these problems by 
using a governing equation in a single parameter, the radial stress. The 
present method is demonstrated to have good convergence characteristics even 
when the initial cr r values are unrealistically large or small. Therefore, 
the nonlinear governing equation in terms of the radial stress and the Newton- 
Raphson technique appear to be the ideal choice for large deflection problems 
of arbitrarily loaded membranes. 


CONCLUSIONS 

Circular membranes with fixed peripheral edges, subjected to arbitrary 
axisymmetric loading, were analyzed. Earlier finite difference methods used 
transverse and radial displacements as parameters and developed two nonlinear 
simultaneous equations. These differential equations were solved simulta- 
neously, to obtain solutions for transverse and radial displacements. However 
because of solution difficulties, only limited results are available. In the 
present formulation a different apprach was taken. Instead of two nonlinear 
simultaneous differential equations, a single nonlinear differential equation 
in terms of radial stress was developed. This nonlinear equation was solved 
by using the finite difference method in conjunction with Newton-Raphson 
method. Three different loading cases, namely, (a) uniformly loaded membrane, 
(b) a membrane with uniform load over an inner portion, and (c) a membrane 
with ring load, were analyzed. 

The results from the present method show that for arbitrary loaded mem- 
brane, transverse and radial displacements differ by large amounts, which 
probably explains why the simultaneous solution method fails to yield a 
satisfactory solution. 
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The present method shows good convergence characteristics for various 
types of loadings* Also central displacement and the central and edge radial 
stresses for uniformly loaded membrane by the present method agreed very well 
with the classical solution. 
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Table l. Sensitivity of the present method 
of the Initial values of a r a 




Number of iterations required 

for convergence 


Uniform loading over the region 

Initial 




c_ values” 





0 < £ < 1 

0 < - < 0.1 

0.5 < - < 0.6 


a 

a 

a 

0.0001 

25 

22 

23 

0.001 

20 

16 

17 

0.01 

14 

10 

11 

0.1 

8 

8 

5 

1.0 

9 

19 

15 

10.0 

20 

27 

26 

100.0 

31 

41 

37 

Converged 




(r - 0) 

0.4310 

0.1104 

0.1237 

Converged 




\ (r - a) 

0.3329 

0.0235 

0.0971 


Membrane idealized with 60 regions. 

^Constant radial stress was assumed at all 61 nodes < 
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TABLE 2. Comparison of Normalized Displacements and Stresses 

for Uniformly Loaded Membrane 


w 





Central transverse 
^deflection 
w at r ■ 0 

Central radial 

stresss 

a at r ■ 0 
r 

Edge radial 
stress 

a at r =» a 
r 

Hencky* 
(ref. 3) 

0.6536 

0.4310 

0.3280 

Kao and Perrone 
(ref. 6) 

0.6541 

0.4289 

0.3306 

Present results 

0.6534 

0.4310 

0.3329 


♦Values taken from reference 6. 
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(a) Circular Membrane 



(b) Loading and deformations 


Figure 1.- Membrane configuration. 
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(a) Uniform loading 



Figure 2,- Types of loading on the membrane. 
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Relative error 





40 


60 8! 


Number of regions, m 


(b) Uniform loading over the region 0 < r/a < 0.5, 
Figure 3 - Convergence study for the membrane. 
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ace merit 







r/a 

Figure 6.- Normalized displacements in a membrane uniformly loaded over 
the inner portion 0 < r/a < 0,1, 

27 
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